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Abstract 

Consider a measure fix = Ylx ^^^^ where the sum is over points x of a Poisson 
point process of intensity A on a bounded region in d-space, and is a functional 
determined by the Poisson points near to x, i.e. satisfying an exponential stabi- 
hzation condition, along with a moments condition (examples include statistics for 
proximity graphs, germ-grain models and random sequential deposition models). 
A known general result says the /UA-measures (suitably scaled and centred) of dis- 
joint sets in M'^ are asymptotically independent normals as A ^ oo; here we give 
an 0(A~^/*-^'^^^-*) bound on the rate of convergence. We illustrate our result with 
an explicit multivariate central limit theorem for the nearest-neighbour graph on 
Poisson points on a finite collection of disjoint intervals. 
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1 Introduction 



There has been considerable recent interest in providing central hmit theorems (CLTs) 
for certain functionals in geometric probability defined on spatial Poisson point processes. 
Such functionals include those associated with random spatial graphs such as the minimal- 
length spanning tree or the nearest-neighbour graph, as well as with germ-grain models 
and random sequential packing models. These functionals are random variables given by 
sums of contributions from points of a Poisson point process in W^. 

A natural extension to random measures may be provided by keeping track of the 
location of each contribution in W^. In this way one can obtain a random field indexed 
by test functions on M"' or by subsets of M.'^. For example, one can consider the measure 
induced by a Poisson process with a point mass at each Poisson point equal to the distance 
to its nearest-neighbour; then a typical multivariate statistic induced by this measure is 
the vector of total edge-lengths of the nearest-neighbour graph on Poisson points over a 
finite collection of disjoint subsets of M^. 

Under certain conditions, it is known [4,10,11] that the measures, appropriately scaled 
and centred, of disjoint sets (or of test functions with disjoint supports) are asymptotically 
distributed as indpendent normals in the large-intensity limit. The object of the present 
paper is to give bounds on rate of convergence; these bounds are the main contribution of 
the present paper. We illustrate our result with an application to the nearest-neighbour 
situation mentioned above. 

The unifying concept of stabilization on Poisson points has proved a useful notion of 
local dependence in the context of geometric probability. This says, roughly speaking, 
that the contribution from a Poisson point is unaffected by changes to the configuration 
of Poisson points beyond a certain (random) distance. 

The methodology of stabilization has been fruitfully employed, in various guises, to 
produce univariate CLTs and laws of large numbers for random quantities in many prob- 
lems in geometric probability; see e.g. [4,8-11,13-16]. The techniques used in this context 
include a martingale method (see for instance [8], and [13] where the method is presented 
for general stabilizing functionals in geometric probability), the method of moments [4], 
and Stein's method [16], which we employ in the present paper. 

The multivariate case, in which several collections of random variables are consid- 
ered, has also received some attention. Applications in geometric probability include, for 
example, the joint normality of certain random spatial graph functionals defined over a 
finite collection of disjoint regions in W^. There are potential applications to multivariate 
statistics, including nonparametric multi-sample tests (see e.g. [17]). 

In the present paper, we employ a form of Stein's method (see [18]), which has the 
advantage that it can provide rates of convergence in the CLT. In this context, Stein's 
method is a useful tool for estabhshing normal approximations and CLTs for sums of 
weakly dependent random variables. In this paper, the weak dependency structure is 
provided by the concept of stabilization on Poisson points. 

In the univariate case, the method yields normal approximation of the sum of a sin- 
gle collection of random variables that are 'mostly independent', i.e. exhibiting a local 
dependency structure. This structure may be captured using dependency graphs. This 
method was first used in the context of geometric probability by Avram and Bertsimas 
in [2] (using the normal approximation error bounds of [3]) to provide CLTs for certain 
random combinatorial structures that are locally determined in some sense, including the 
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j'-th nearest-neighbour graph, and the Delaunay and Voronoi graphs. 

Using the sharper normal approximation bounds of [5], more general results for uni- 
variate normal approximation based on Stein's method for random point measures were 
given by Penrose and Yukich in [16]. That paper is the foundation for the present work, 
which is its multivariate analogue. 

Multivariate CLTs for random measures in geometric probability have recently been 
proved via the method of moments [4] and also the martingale method [10]. In particular, 
[10] also covers lattice processes (such as percolation), and does not require 'exponential' 
stabilization, and so admits a larger class of measures. The advantage of the results in 
the present paper is that information on rates of convergence is provided. 

Beyond the context of geometric probability, mulivariate central limit theory has been 
well studied. Related results include multivariate central limit theorems for sums of in- 
dependent random variables given in [6]. In [7,17], multivariate normal approximation 
bounds are given for sums of (locally) dependent random variables, often chosen in some- 
what special ways, including certain statistics defined on random graphs. The results 
in the present paper have the advantage of being more generally applicable in geometric 
probability. 

2 Main result 

The basic setting follows that of [16]. Let d e N. As in [16], we consider marked point 
processes in M.^ for the sake of generality. Let {A4,J^M:^m) be a probability space (the 

mark space). Let ^(.x, s; X) be a measurable [0, cxD)-valucd function defined for all triples 
{x, s; X), where x G M'', s E Ai arc such that {x, s) G X, where X G M."^ x Ai is finite. 
When {x, s) G (M"^ x A4) \ X, we abbreviate notation and write ^{x, s; X) instead of 
C{x,s;XU{{x,s)}). 

Given X C R"^ x M, a > and y G M!^, set y + aX := {{y + ax, s) : {x, s) G X}, 
i.e. translation and scaling act only on the 'spatial' part of X. For all A > let 

^x{x, s; X) := ^{x, s-x + A^/''(-x + X)). 

Thus is a 'scaled- up' version of ^, defined on a scaled- up version of the (marked) point 
set X dilated around x. We say that ^ is translation invariant if C,{x + y, s;y + X) = 
^{x,s;X) for aU y G M'^, aU {x,s) e M."^ x M and all finite X C M."^ x M. When ^ is 
translation invariant, the functional simplifies to ^x{x, s; X) = ^{X^^^x, s; X^^'^'X). 

For q G [1, oo], let || • \\q denote the £q norm on M.^. In the sequel we will use q — 2 
(the Euclidean norm) and q = oo. 

Let K be a probability density function on M'^ with compact support A C M'^, where 
A is non-null (i.e. has non-zero Lebesgue measure). We assume throughout that k, is 
bounded with supremum denoted by ||«;||oo < oo- For all A > let V\ denote a Poisson 
point process in R'^ x Al with intensity measure {XK{x)dx) x Pjv/((ds). 

Wc use the following notion of exponential stabilization, as given in [16] (taking the 
A\ there to be A for all A). For x E M.'^ and r > 0, let Br{x) denote the Euclidean ball 
centred at x of radius r. Let U denote a random element of A4 with distribution P^k, 
independent of V\. 
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Definition 2.1 ^ is exponentially stabilizing with respect to k and A if for all X > 1 and 

all X E A, there exists an almost surely finite random variable R := R{x, A), (a radius of 
stabilization for ^ at x) such that 

^(x, u- [Vx n {b^-,/,r{x) X M)] yjx) = u- Vx n {b^-,„r{x) x m)), 

for all finite X d {A \ By^-i/dji{x)) x ^A, and moreover 

limsup t'^ log ( sup P[i?(x, A) > t] ) < 0. 

Roughly speaking, A) is a radius of stabilization if the value of ^x at x is unaf- 
fected by changes to the configuration of Poisson points outside By^-i/d^{x). Exponential 
stabilization is known to hold for many 'locally determined' functionals defined on spatial 
point processes, and in particular in several cases of interest in geometric probability; see 
for example [16]. Following [16], we also make the following definition. 

Definition 2.2 ^ has a moment of order p > (with respect to k and A) if 

sup E[|eA(x,f/;PA)n<oo. (2.1) 

x>i,xeA 

For A > 0, we define the random weighted point measure /i^ on R"^, induced by ^a, by 

:= ^xix,s;Vx)6x, 

{x,s)(^Vxn{AxM) 

where 5^ is the point measure at x G M'^. 

For r C M.'^, let B{r) denote the set of bounded Borel-measurable functions on F. For 
/ G B{T), let (/, /i^) := /p fdfiy Let $ denote, as usual, the standard normal distribution 
function on R. We recall the following univariate normal approximation result of Penrose 
and Yukich (contained in Theorem 2.1 of [16]). 

Proposition 2.1 [16] Let ^ be exponentially stabilizing and satisfy the moment condition 
^2. 1\) for some p > 3. For F a non-null Borel subset of A, let f G B{T) and put T : = 
(/, Then there exists a constant C G (0, oo) depending on d, ^, f, and k such that 
for all \>2, 

<C(logA)3^A(Var[r])-='/2_ ^2.2) 

For fixed m G N, let Fj, z = 1, ... ,m be non-null Borel subsets of A C M'^. For 
notational simplicity, for i = 1, . . . ,m and for /j G B(Ti) set Tj := {fi,^\) = fid^y 
These are the quantities of interest to us in the present paper. By Proposition I2.ll under 
appropriate conditions, we have that, individually, each Tj satisfies a normal approxima- 
tion result of the form of (12.21) . For the present paper, we will impose one extra condition 
to control variances such as Var[Tj]. 

(Al) There exist constants Cj G (0, oo) such that for each i, for all A sufficiently large, 
Var[T,] > QX. 



sup 



P 



■ T - E[T] 
(Var[T])i/2 



< t 
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Under assumption (Al), the bound on the rate of convergence on the right of fl2.2p (in the 
case T = Tj) becomes 0(A~^''^(log A)^'^) (compare Corollary 2.1 of [16]), and in particular 
(12. 2p yields the central limit theorems 



(Var[T,])i/2 



V 



Ar(0,l), 



as A — ^ oo, where A/'(0, 1) is the standard normal distribution on M and ' — denotes 
convergence in distribution. Condition (Al) is true in many cases. In Section WA\ we will 
give some sufficient conditions for (Al) to hold, and discuss alternative conditions which 
lead to somewhat stronger versions of (Al). In particular, it is often possible to show 
(under appropriate conditions) that A~^Var[Tj] af for some af G (0, oo), which may be 
'explicit' (see Section B?T1) . 

Our main result. Theorem 12.11 below, extends Proposition 12.11 to give a multivariate 
central limit theorem for (Tj : i = 1, . . . ,m), centred and scaled, with a bound on the 
rate of convergence. We impose the additional assumptions that (Al) holds and that 
the sub-regions Fj are pairwise disjoint and satisfy the natural regularity condition (A2) 
below. The central difficulty in extending Proposition 12 . 1 1 to a multivariate version is that 
the Tj are not, in general, independent. However, with the aid of stabilization we will 
show that they are 'asymptotically independent' in an appropriate sense. 

To state (A2), we introduce some notation. For measurable B C M*^, let \B\ denote the 
((i-dimensional) Lebesgue measure of B. Let dB denote the boundary of B. For B cM.'^ 
and X eR'^ let dq{x,B) := inf \\x - y\\q. Also, for B,B' C M"' with 5 n 5' = 0, let 
dq{B,B') := misceB,yeB' \\x — y\\q, i.e. the shortest distance (in the iq sense) between B 
and B'. For r > 0, let dr{B) denote the r-neighbourhood of the boundary of i? C M'^ in 
the norm, that is the set {x G M'^ : d^{x, OB) < r}. 

(A2) For each \dr{Ti) \ = 0{r) as r j 0. 

Sufficient conditions for (A2) include that each of the Fj is convex, or each is the finite 
union of convex regions (e.g. polyhedral). We can now state our main result. 

Theorem 2.1 Let C, be exponentially stabilizing and satisfy the moment condition Ii2. 1\) 
for all p > 1. Let m G N. Let Fi, F2, . . . , Fm be fixed disjoint non-null Borel subsets of A 
satisfying (A2). Fori = 1, . . . ,m, let fi & 'S(Fi) and setTi := {fi,fi^^). Suppose that (Al) 
holds. Let e > 0. Then there exists a constant C G (0, 00) depending on d, ^, k, e, {/j} 
and {Fj}, such that, for all A > 1, 



In particular, from (12.31) we obtain the multivariate central limit theorem that says 



as A ^ 00, where A/'(0, J^) is the m-dimensional normal distribution with mean and 
covariance matrix given by the identity matrix J^. It was already known [10,11] that under 




(2.3) 




(2.4) 
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similar conditions to those of Theorem 12.11 we have fl2.4p . at least when A~^Var[Tj] af 
for some af G (0, oo). Theorem 12.11 adds to this by providing a bound on the rate of 
convergence. 

As an example of the application of Theorem 12.11 one can take fi = Ir^ for i = 
1,2, ... ,m, where Ir is the indicator function of F C M'^. We indicate some particular ap- 
plications of Theorem 12. II in Section HI Under additional technical conditions, one can say 
more about the asymptotic behaviour of the variance terms in (12.31) : see Section |4?T] below. 

Remark. The relatively slow rate of convergence in higher dimensions arises primarily 
due to the possibility of strongly dependent points in the neighbourhood of the interface 
of adjacent regions. If all of the Fj are separated by a strictly positive distance, then 
our methods can be adapted to yield a rate of convergence of the same order as in the 
univariate result (Proposition 12. ip . that is 0(A~^/^(log A)^"^). 

For ease of presentation, we prove Theorem 12. 11 in Section [3] under the conditions that 
^ is translation invariant and that the mark space is degenerate (i.e. A4 = {1}), and so 
from now on we suppress any mention of A^. In particular, point sets such as X and V\ 
will be treated as (their corresponding) subsets of M"^, and we will write C,x{x', P^) rather 
than ^\{x, 1; X). The proof can be adapted for the general marked case, as in [16]. 

3 Towards a proof of Theorem 12.11 

For everything that follows, we assume that Fi, F2, . . . , F^ are (arbitrary) non-null Borel 
subsets of the bounded region A C W^, such that Fj fl Fj = for i ^ j, and condition 
(A2) holds. Also, for each i we have a function fi e B{Vi). 

For fixed a > 0, let sx := aA'^/'^log A, and let F^'^ denote the sx 'boundary region' of 
Tj ^ A, in the sense 

Tf:= [x eV,: d^{x,dV{) <aX-^/Hog\] = V,f\d,^{Vi). (3.1) 

The remainder of the set Fj we simply call the 'interior' and denote by F™, where 

Vf ■= [xeVi- d^{x,dVi) > aA-^/'^logA} = Fi\5s,(Fj). 

As previously mentioned, we assume that ^ is translation invariant, and that Ai = {1}. 
Define 

T^^ := [ fdfih and := / /.d/,, 

I % 

SO that Tj = + Tf^. To prepare for the proof of Theorem 12. II we need some auxiliary 
lemmas. For the subsequent results, we will need the following covering of scaled-up Borel 
regions X^^'^B C by cubes of side 1. 

First we need some more notation. Let card(A') denote the cardinality of set X. For 
X G M*^, let Qx denote the unit- volume ioo ball in with centre x (i.e., the unit d-cuhe 
at x). For a Borel set B C A C R"^, let 

ZxiB) := {xeZ'^-.Q^n X^/'^B ^ 0} , (3.2) 
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and set nx{B) := caid{Zx{B)). Then the covering of X^^^B is 

Qx{B) := {Q,:ze Zx{B)}. (3.3) 

The next resuh gives error bounds for approximating the volume of X^^'^Ti or of A^^'^F^'^ 
(as defined at (13.11) ) by the number of unit cubes in Z"^ in its covering (as defined at (13.21) 
and dSJD). 

Lemma 3.1 Let Fj be a non-null Borel subset of A G such that \dr{Ti)\ = 0(r) as 
r I 0. Then, as X ^ oo, 

nx{T,)-\X'/%\=0{X'^'-'y'). (3.4) 

Define F^*^ as at ^3. 1]) . Then, as A — > oo, 

nx{Tf) - I A^/'^F^'^I = O (A^'^-^)/'' log A) . (3.5) 

Proof. There exists a constant c G (0, oo) (depending only on d) such that, for any 
A > 0, and any non-null Borel subset B of A, 

X^/dQ c y C X^''^B U d^{X^''^B), 

z(iZ^{B) 

and hence 

\X^'''B\ < nx{B) < \X^/''B\ + \d,iX'^''B)\ = \X'^''B\ + X\d,x-i/4B)\. (3.6) 

In the case B = Fj, the regularity assumption that |(9r(Fj)| = 0(r) as r | implies that 
l-^cA-v^r^)! = 0(A-i/'^). Thus dilD follows from flgTBjl . 
In the case B = F^*^, we have that 

|9,,-v.(rr)l<|9,,-v.+.,(r.)l = o(s,), 

as A — ^ oo, again by the regularity assumption on Fj. Thus (13. 6p yields (13. 5p in this case. □ 

Once more consider a Borel subset B of A C M.'^ and the covering Qx{B) of X'/'^B. 
For all z G Zx{B), the number of points of Vx H X^^^^Qz is a Poisson random variable 
Nz with parameter z/^ := A Jx-i/dq K,{x)dx. Assuming Uz > 0, choose an ordering on 
the points of Vx H X"^^'^Qz uniformly at random from all Nz\ possibilities. List the 
points as i, . . . , X^, ^r^, where conditional on the value of Nz, the random variables 
Xz,k, k = 1,2, Nz are i.i.d. on X^^^'^Qz with a density Kz{-) ■= k{-)/ Jx-i/dg K{x)dx. 
Thus we have the representation 

VxnB= U [j{{Xz,k}nB). 

z€Zx(B) k=l 

Then for / in B{B), we can express (/, /x^) as follows: 

zeZxiB) k=l 
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For all z G Z\{B) and for all A; G N, let Rz^k denote the radius of stabilization of ^ at 
^z,k if 1 < < and let Rz,k = otherwise. Define the event E^^k '■= {Rz,k < a log A}. 
We define here the function T{B; f) as follows, the idea being that T{B] f) is, with high 
probability, the same as (/, /x^), but exhibits a much more localized dependency structure. 
Set 

f{B-f):= J2 (3.8) 

zeZxiB) k=l 

where we use to denote the indicator random variable of the event E. 

Recall that Fj, i = l,2,...,m are disjoint non-null Borel regions in A C M'* and 
fi G B(Ti) for i = 1, 2, . . . , m. Then for each i, T(Fj; /,) is defined by (13.81) . In the same 
way as we use the abbreviations Tj, T}"^ and T™, we let Tj := T(Fj; fi), T}"^ := T(F^^; fi), 
and fi'' := T(F^"; fi). Thus fi = f^'^ + f^. 

For z G Z\{B) let Yz{B] f) be the contribution to T{B] f) from the points in X'^^'^Qz, 

i.e. 

Yz{B-f) := Y,UXz,k\Vx) ■ 1e^,, ■ f{Xz,k) ■ 1b(X,,,), (3-9) 
fc=i 

so that f{B;f) = j:^^^^^^^Yz{B;f). 

Let Aa, A > 1 be a family of Borel subsets of A C M*^. The next two results show that 
the moments condition (12. ip implies bounds on the moments of Yz{A\; f) for / G B{A). 
When we come to apply the two lemmas below, we will be taking Ax = Fj or Ax = T]"^. 

Lemma 3.2 Let Ax, X > 1 be a family of Borel subsets o/ A C M°'. // Ii2.1\) holds for 
some p > 0, then there is a constant C G (0, oo) such that for all A > 1, all k > 1 and 
z G Zx{Ax) 

Y\\ix{Xz,k\Vx) ■ lA,(X,,fc) • l{fe<jv4r] < C. (3.10) 

Proof. It suffices to consider the case with Ax = A for all A. The proof of the 
lemma closely follows that of Lemma 4.2 in [16], although our covering is somewhat 
different. In the notation of the proof of Lemma 4.2 of [16], we have pa = 1 and 
Vi = yz, = X fx-i/dQ^_ K(x)dx < ||k||oo, where we have written Zx{B) = {zi, . . .,Zn^^B)}- 
Then, following the argument in [16], we obtain (13.101) . □ 

Lemma 3.3 Let Ax, X > 1, be a sequence of Borel subsets of A G M°', and suppose 
f G B{A). If /12. 1\) holds for some p > 1, then for any q G (l,p) there is a constant 
C G (0, oo) such that for all X > 1 and all z G Zx{Ax) 

\\YziAx;f)\\l<C. (3.11) 

Proof. The proof closely follows that of Lemma 4.3 in [16], again with px there equal 
to 1 (and Ui < ||k||oo)- Thus, with the use of Lemma [3^ (and the boundedness of /), we 
obtain (ISH]) . □ 
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Lemma 3.4 Suppose that ^ is exponentially stabilizing and satisfies the moments con- 
dition Ii2.1\) for some p > 3. Then there exists a constant C G (0, oo) such that for all 
A > 2 

sup 

Moreover, Ii3.12\) holds with % replaced by T-^ everywhere. 

Proof. The statement for Tj follows from equation (4.18) in [16] with px = O(logA), 
g = 3, and taking the Ax of [16] to be Fj. In equation (4.18) of [16], T( is the equivalent of 
our Tj, Tx is our Tj, and S is our (Tj — E[Tj])(Var[Tj])^-'^/^. The statement for follows 
in the same way, this time taking the Ax of [16] to be Tf^. □ 

Lemma 3.5 Suppose that h2. 1\) holds for somep > 2. Then there exist constants Ci, C2, C3 G 
(0, 00) such that, for all A > 2, 

Var[f^'^] < CiA('^-i)/'^(logA)'^+\ (3.13) 
Var[fi] < C2A(logA)^ and (3.14) 
Var[T;'°] < C3A(logA)'^. (3.15) 

Proof. First we prove fl3.13p . Consider the covering QxiT'^'^) of X^/'^T}'^ by unit rf-cubes, 
as defined at (fSJ]). For z G ^A(r^'^) let Y.iT^'^; fi) be the contribution to f^'^ from the 
points in X'^^'^Qz, as defined at (13.91) . that is 

y,(F^^; /,) := J2UXz,k;Vx) ■ 1e^,, ■ /.(X,,^) ■ lrbd(X,,fc). (3.16) 

k=l 

Now, using the representation T}"^ = J2zez^{r^'^) '^ziX}'^', fi), we have 

V8.T[f^'] = ^Var[F,(Fr;/,)] + 5^Cov[F,(Fr;/,),FUFr;/,)]. (3.17) 

z z^w 

By the assumption that (12. ip holds for some p > 2, by taking q = 2 and Ax = F^*^ 
in Lemma 13.31 we have that Var[Y2(F^'^; /j)] < V, for some constant V < 00, for all 
z G Zx(T]"^). So by the Cauchy-Schwarz inequality we have Cov[F2(F^'^; /«), Yy^iT^'^; fi)] < 
V. Also, Yz{T]"^;fi) and F^(F^^;/i) are independent if d2{Qz,Qw) > 2a log A (by the 
definition of Ez^k)- Further, given z, the number of w for which d2{Qz,Qw) < 2a log A is 
0((logA)'^). Hence dXTT]) implies that 

Var[T;^'i] < nx{T^'){V + 0{{logX)'')). (3.18) 

Then by (13. 5p we have that 

nxiT'^'') = AIF^-^I + O (A('^-i)/'^ log A) = 0{X^''-'^/'' log A), (3.19) 

using ([31]) and (A2). So from fl3l8D and fl3l9|) we obtain fl3J3D . 

The proof of (13.140 follows similarly, using Ax = Fj for all A in Lemma 13.31 and (13.40 
in place of (13.51) . Finally, (13.151) follows from (13.141) . (I3.13P and the Cauchy-Schwarz 
inequality, since T™ = % — T}"^. □ 
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P 



- E[f,] 
(Var[f,])i/2 



< t 



<CA(Var[T,])-3/2(logA) 



\3d 



(3.12) 



Lemma 3.6 Suppose that ^ is exponentially stabilizing and satisfies the moments condi- 
tion Ii2.1\) for some p > 3. Then there exists a constant C G (0, cxo) such that for any 
(5 > 0, all X> 2, and any t eM. 



P 



(Var[r,])i/2 



< 6 



< \l =-6 + c {log \f'\(y^T[f,]r'/\ 



(3.20) 



and also 



P 



(Var[T,])V2 



< 6 



< 2^/^(5 + C(logA)3'^A(Var[fi])-3/2 



+P 



J-bd _ ^rp^ 



bdl 



(Var[r,])i/2 

Proof. First we prove fl3.20p . For the duration of this proof, write 



> 6 



(3.21) 



F{t) = P 

Then we have that for t G M and 6 > 



(Var[r,])i/2 



< t 



P 



t 



< 6 



F{t + 5)- F{t - 5) 



(Var[r,])V2 

= $(t + 5)- ^t -5) + [Fit + 5)- ^t + 5)] - [Fit -5)- $(t - 5)] 
< |$(t + 6)-^t-6)\ + \F{t + 6)-^t + 6)\ + \F{t - 5) - ^(t - 5)| . 



Then (13.201) follows from the Mean Value Theorem (applied to the first term on the right 
of the above inequality) and Lemma 13.41 (applied to the other two terms). Finally, we 
have that for 6 > 



P 



fin _ ^fii 



< 6 



< P 



(Var[T,])V2 
Then using (^M> yields (Km . □ 



(Var[T,])V2 



< 26 



P 



bdl 



(Var[T,])V2 



> 6 



Lemma 3.7 Suppose that the moments condition ^2. ij) holds for allp > 1, and condition 
(A2) holds. Let k be an even positive integer. Then there exists a constant C G (0, cxd) 
(depending on k) such that for all X> 2, 



E 



^bd _ ^j.^ 



bdl 



< (j^k{d-l)/{2d)(j^^^^kil+d)/2_ 



(3.22) 



Proof. Again consider the covering Qx{T^'^) of A^/'^F^'^ as defined at (13.30 . For z G 
ZA(r^^), let be the contribution to f}"^-K[f}"^] from cube Q„ that is := Y,{T^'^; fi)- 
E[Y,{T^^-Ji)] where Y,{T^'^; fi) is given by (lAj) . Thus, for all z G ZxiT^"^), E[Y,] = 
and Var[y2] < V for constant V, by Lemma [3.31 
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Let k be an even positive integer. Then 

The term E [Yz-^Y^^ ■ ■ ■ ^z^] will vanish if any of the cubes corresponding to the Y^j is 
farther than 2a log A from all the other cubes (since then it will be independent of the 
other Yz- and has expectation zero). In other words, the term vanishes if the appropriate 
geometric graph (in the sense of [9]) on zi, Z2, . . . , Zk has any isolated vertices. For a 
non-zero contribution to the sum, we require the graph to have no isolated vertices — 
so it must have no more than k/2 components. So in effect, there are at most k/2 'free' 
indices of {zi, . . . , Zk). Values that are not 'free' have 0((log A)*^) possible values. 

Further, E [F^^F^j ' ' ' < C for some constant C, by Lemma [3^ (given the moments 
condition (12. ip for all p > 1) and Holder's inequality. Thus for some other constant also 
denoted C, 

E E E ^[YzX.---Yz,] < C{n,{Tf)r/'{\ogXr/' 

the final inequality by (13.51) . (13.11) and (A2). Hence we have (13.221) . □ 



E 



T 



bd 



ET 



bd 



The next lemma says that given condition (Al), we can obtain lower bounds on the 
variances of and Tj. We will need the following result from [16] (see (4.17) therein), 
which says that if ^ is exponentially stabilizing and satisfies the moments condition (12.11) 
for some p > 2, then 



Var[f,]-Var[r,] 



< CA" 



(3.23) 



Lemma 3.8 Suppose that (Al) and (A2) are satisfied, and that the moments condition 
Ii2.1\) holds for all p > 1. Then there exist constants C G (0, cxd) and Aq G [1, cxd) such 
that for all A > Aq 



Var[T,] > CX, 
Var[T;'"] > CX. 



and 



(3.24) 
(3.25) 



Proof. These follow in a straightforward manner from (I3.23p . (Al), (I3.13p . (I3.14p and 
the Cauchy-Schwarz inequality. □ 

Lemma 3.9 Suppose that ^ is exponentially stabilizing and satisfies the moments condi- 
tion \2.1\) for all p > 1. Suppose conditions (Al) and (A2) hold. Then for any e > 0, 
there exists C G (0, oo) such that for all X > 1 



< CA-^/(2d+.)^ 




- E[T,] 
(Var[Ti])i/2 

fin _ ^fin 

(Var[T,])i/2 



1=1 

m 



i=l 



(3.26) 
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Proof. We abbreviate our notation for the duration of the current proof by setting 
(Tj := (Var[Tj])^/^. Then we have 



P 



< 



P 



i=l J 1=1 

m ~\ m 

f]{(Tr-E[rr\}ar^<t^ -Hhu) 

i=l J 1=1 

+ ^P \(fr - E[rr\)a-' < U, (T, - E[r,])ari > U 

i=l 
m 



(3.27) 



1=1 



For any /3 > 0, we have 



max f P 



p 



< p 



>x-(' 



p 



iTr-E[Tn)a-'-t 



. (3.28) 



Then, from (KTlh and fl3:28|) 



P 



< 



P 



f|{(T.-E[T,])ari<t, 

i=l 

[]{{fr-E[Tr\)^r^<u 



i=l 
m 



i=l 



+2^P [|(7;^^-E[f;'^<^])(Tri >A~'^j +2^p[|(f;'°-E[f;'°])ari-tJ < A-'^j .(3.29) 



i=l 



i=l 



Since (i2(A^/'^r™, A^/'^F™) is at least 2a log A for i ^ j, T™, 1 < i < m is a sequence of 
mutually independent random variables, so that 



P 



m 

f]{ifr-E\fr])ar^<U 

i=l 



l[F\(fr-E\fr\}crr^<U 



(3.30) 



1=1 



Also, from Markov's inequality, we have that, for k G 2N, 



P 



iT^''-E[T,^^])a, 



> 



A-^ 



< E 



bdi 



~ \ -k/2 

Var[T,]j A^^. (3.31) 



Then we obtain, from (ICTj) . with (Km and (Km . 



P 



(T^^-EiTrM 



bdi\„-l 



> 



A-^ 



^^^kW-l/i2d))^^^^^^k{l+d)/2. 



(3.32) 
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this then gives a bound for the penultimate sum in fl3.29p . To bound the final sum in 
(^M), taking 6 = A"^ we have from ([MID, fl332|) and iKm that 



P 



{Tr-E[Tn)ar'-t 

< 2^X-^ + C(l0g A)3'^A-^/2 ^ ^;^fc(/3-l/(2d))(i^g;^)fc(l+d)/2 



(3.33) 



To obtain the best rates of convergence via this method, we want to maximize the lowest 
power of A~^ on the right-hand sides of (13.321) and (13.331) . So we choose j3 such that 
-P = k{p- l/(2c?)), that is, take 



k 



2d{k + l) 



(3.34) 



For any £ > we can choose k large enough in (13.341) to give l/{2d) >/?>!/ {2d + e/2). 
Then, for A sufficiently large, A'^/^^d+e) > ^-i/{2d+e/2)QQg^)fc(i+d)/2_ ^^^^ f|3:29D and 

( 13.301) . with the bounds (13.321) and (13.331) we obtain (I3.26p . This completes the proof of 
the lemma. □ 

Proof of Theorem 12.11 To complete the proof we proceed in a similar manner to [16]. 
Let 

m Nz 

i=i zeZx{r^) k=i 

recalling the definition of the event Ez^k just below (13.71) . By standard Palm theory 
(e.g. Theorem 1.6 in [9]) and exponential stabilization (see (4.11) in [16]), we have that 
IP[-E'^] < CX"^ for a sufficiently large and some C G (0, oo). Then |Tj — Tj| = except 
possibly on the set E^, which has probability less than CA^^. 

For ^ = 1, . . . ,m, let := (Var[T,])-i/2(T, -E[T,]) and Z, := (Var[T,])-i/2(T, -E[T,]). 
Then for 5 > we have that for any e M 

{{Z, < U)A{K, < U)} C{\K,- U\ <6}U{\Z,- K,\ > 5}, 

so that 



P 



n {Z^ < ti} 



.1=1 



< 



p 



i=l 



i=l 



+ Y,^[\K,-t,\<6] + Y,^[\Z,-K,\>6] 



(3.35) 



i=l 



i=l 



Then, using ( I3.26P for the first term on the right-hand side of the inequality in (I3.35p . 
and (I3.20p with fl3.24p for the second, we obtain 



P 



i=l 



i=l 



< CA-^/(2d+e) ^ 



i=l 



< u 



i=l 



(Var[T,])i/2 

m 

+C5 + C(log XfX-^'^ + ^ P [|Z, - ir, I > (5f3.36) 



i=l 
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We now consider the second term on the right-hand side of fl3.36p . For ease of notation, 
write 



(Var[f,])V2 



< t 



for z = 1, . . . , m. Then 

m rn 



i=l 



i=l 



i-1 



j=i+l k=l 



i=l 



(3.37) 



1=1 



Writing 



H,{t) :=P 



fm _ ^fh 

(Va^[T^/2 



< t 



we have that, for i = 1, . . . , m 



\G^{U) - < \H,{Uil + 7,)) - + 7,))| + |$(t,(l + 7,)) - , (3.38) 

where 1 + 7^ := (^ y^^jfil-^ ^ ■ Then, using Lemma [13] we have that the first term on the 
right-hand side of (13.381) satisfies 

mti{l + ^,)) - $(t,(l + 7,))| < C(logA)3'^A(Var[f;-])~3/2 < C\-'/\\og\)''', (3.39) 

by (13.251) . In order to deal with the second term on the right-hand side of (13.381) . we need 
to estimate •ji. We note that 



VarJ^L 
Var[f;''^] 



1 + 



Var[7;i°] Var[7;'^ 



Then using the upper and lower variance bounds (13.131) . (I3.15p . (13.251) . and the Cauchy- 
Schwarz inequality, yields 

= 1 + 0(A-i/(2^)(log A)(2'^+i)/2) 



so that 



7, = 0(A-l/(2ci)(lQg^)(2d+l)/2) 



(2d+l)/2\ 



(3.40) 



Since for all s < t we have |$(s) — $(t)| < {t — s) sup^<u<( vi.'^) (where (p is the standard 
normal density function), we have 

sup|$(ti(l + 7i)) -$(^01 



< Csup |t,|A-i/(2'^)(logA)(2^+i)/2 

ti \ 

< CA-l/(2d)^i^g;^)(2d+l)/2_ 



sup V9(m) 

|«-iil<iiCA-l/(2d) (log A)(2d+l)/2 



(3.41) 
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So, for the second term on the right-hand side in fl3.36p . we obtain from fl3.37p . fl3.38p . 
( Km and K^ 



sup 

ti,...,t„ 



np 

i=l 



fin _ ^fin 

(Var[f,])i/2 



< u 



i=l 



(3.42) 



< C(l0gA)=^'^A-l/2 ^ ^^-l/(2d)(i^g^)(2d+l)/2^ 

We now move on to the fifth term on the right-hand side of fl3.36p . We have 

\Z,-K,\ = (Var[T,])"i/2|(T,-E[r,])-(T,-E[T,])| < {Va.i[f (^|t, - r,| + E[|T, - T,|]) , 

and from just below (4.19) in [16], we have that this is bounded by CX^^ except possibly 
on the set which has probability less than CX~^. Thus by (13.361) with 6 = CX~^, and 
using (I3.42P for the second term on the right-hand side of (13.361) . we obtain 



sup 

ti,...,t„ 



P 



.i=l 



i=l 



< CA-i/('''+^) +C(logA)3''A 



-1/2 



+CA-l/(2d)(iQg^)(2d+l)/2 ^ ^^-3 ^ o(^-l/(2d+.)) 

By the triangle inequality we have 



(3.43) 



sup 

^1 V ■i^m 



P 



n 



< sup 

tl,...,trr 



P 



T.-E[T.] / Var[TA 
|J\(Var[f,])V2 - ^ \YaT[fJ 



Ml(Var[T,])V: 

1/2 



■1=1 



+ sup 

t\,...,t„ 



i=l 



i=l 



.(3.44) 



Now from (13.230 and (13.241) . there is a constant C G (0, oo) such that for all A > 1 and 
all h e R 



/Varp^A 
VVar[r,]; 



1/2 



It; I 



^ ^ Var[T.]-Var[f,] \ _ ^ 
Var[T,] y 



(l+0(A-3))'/'-l <C|ti|A 



-3. 



then since for all s < t we have |$(s) — $(t)| < (t — s) maxs<u<t we get 
'Var[T,]^^/'' 



sup 



< C sup I |tj|A 

ti 



sup (^(u) < (:7A"'^(3.45) 



,Var[r,], 

Then, considering the second term on the right-hand side of (I3.44p . we have 



sup 

^1 )• • • 



n-('-(Si)")-n» 
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sup 

tl,...,t„ 



< 



sup 



tl,...,t 



i=l 



VVar[T,]; 
VVar[f,]; 



(3.46) 



by (K4E\\ . Thus for any e > 0, from and (13:431) with ([3361), 



sup 

^1 ) ■ ■ ■i^m 



P 



li(Var[T,])i/2 



This completes the proof of Theorem I2.11 □ 



4 Indication of applications 



In applying Theorem 12. 11 one needs to check that the stabilization and moments conditions 
given in Definitions 12.11 and 12.21 hold. These conditions, or related versions thereof, are 
known to hold for many problems of interest in geometric probability; see [4] and [16] 
for an indication of problems for which exponential stabilization and moment bounds are 
satisfied. 

One also needs to verify the variance bound (Al): we discuss methods of doing this 
in Section W7\\ below. In many cases, (Al) (or related versions thereof) has been demon- 
strated, see for example [2,4,13]. 

In Section 14.2! we give an example of our result as applied to the A;- nearest neighbour 
graph. In particular, we give a multivariate CLT with explicit variance scalings in the case 
of the nearest-neighbour (directed) graph on disjoint subsets of the real line (Theorem 14.11 
below). 



4.1 Control of variances 

In this section we discuss conditions under which one can say something about the vari- 
ances Var[Tj]. Recall that Theorem 12. H is stated under assumption (Al). First we give a 
sufficient condition for (Al) to hold, similar in spirit to that used by Avram and Bertsi- 
mas [2]. Once again, for notational convenience we consider only the unmarked case with 
M = {1}. 

First we introduce some notation. Recall that denotes the unit d-cube centred at 
X G W^. For a non-null Borel subset B oi A d and A > we define the following 
packing of X^^'^B by unit d-cubes. For X^/^'-B C let 

Wx{B) := {weZ'':Q^(Z X^'^B] , (4.1) 

and set m\{B) := card(WA(-B)). Then we define the packing lCx{B) of X^^'^B by 

ICxiB) := {Q^:we Wx{B)}. (4.2) 

Let / e B{B). For w G Wx{B) set 

F^:=Fx{Q^;B):= ^ 6(x; Pa) ■ /(x). (4.3) 

Let J^x denote the a- field generated by the points of Vx- 
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Definition 4.1 Let {A^ : w e W\{B)} be a set of mx{B) events in T\, associated with 
the m\{B) cubes Qw{B), so that each A^^ occurs with probability uniformly bounded away 
from zero. Let J = {wi, . . . , wm} be the (random) set of indices w e yVx{B) such that 
occurs. Let Q\ denote the a -field generated by the random set J C >Va(-B) and the values 
ofFx{Q^-B)forwiJ. 

We say that fi^ nondegenerate on {B, f) if there exist events {A^j : w G Wa(-B)} in 
T\, with P(A^) > p > for all w, such that: 

(i) given Q\, for all w & J, Var[F^|^A] > ''7 > a.s.; 

(a) given Qx, for all w,v & J with w ^ v, and Fy are (conditionally) independent. 

The idea of this condition is that the events A^ essentially 'isolate' cubes Q^, while 
allowing strictly positive variability (of the integrated measure) within the cube, and a 
positive fraction of all the cubes Qw will be so 'isolated'. 

This nondegeneracy condition can often be demonstrated. In many cases, event A^j 
will involve a configuration of many points in an 'annulus' just outside the cube Qw, and 
an empty 'moat' inside the cube, that ensures sufficient independence; see [2] for such 
a construction (in a similar context) for the total length of the j-th nearest-neighbour, 
Voronoi, and Delaunay graphs. 

We now show that given the nondegeneracy condition of Definition 14.11 we have lower 
bounds of order A on the variances of Tj. 

Lemma 4.1 Let T be a non-null Borel subset o/ A C M*^ such that \dr{T)\ = 0{r) as 
r I 0. Then, as X ^ oo, 

mxiV) - \\^/'^T\ = O (a(^-^)/^) . (4.4) 

Proof. This follows in a similar way to the proof of (13. 4p given previously. □ 

Lemma 4.2 Suppose that /x^ is nondegenerate on (Fj, /j) (see Definition ^.i| ). Then there 
exists a constant Ci G (0, oo) such that for all A sufficiently large 

Var[Ti] > CiX. 

Proof. From the definitions of the packing and covering defined by (14. ip . (14. 2 p and (I3.2p . 
(13. Sp respectively, and the equations (13. 7p and (14. 3p . we have that for Borel F C A C M*^ 
and / G B{T) 

U.A)= E i^A(Q.;r) + A,(F), (4.5) 

where we have set 

Aa(F) := Y.^x{Xu.,k\Vx)-f{Xu,,k)-lT{X.^,k). 

w&Z),{T)\W,,{V) k=l 

That is, Aa(F) gives the contributions to (/, p^) from cubes that are in the covering of 
X^/'^Y but not the packing. 
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Consider F = Fj with |Fj| > and f = fi G B{Ti). By a similar argument to fl3.13p . 
and fl3.23p . we have that Var[AA(Fj)] = o(A) as A ^ oo. So, by (14.51) and the Cauchy- 
Schwarz inequahty, to prove the lemma it suffices to show that for all A sufficiently large 



Var 



> CX, 



for some C G (0, oo). 

Recall Definition 14.11 The proof now follows the idea of Avram and Bertsimas (see [2], 



Proposition 5). Let M = ^ 



Recall the packing defined by (14. 2p . Index the 



(4.6) 



cubes Qw for which A^, holds by J = {wi, . . . , wm} ^ y^xi^i)- Then 
E[M] = ^ P(A^) > pm^{T,) > CX\T,\ > CX, 

using (14. 4p for the penultimate inequality, and the fact that |Fj| > for the final one. As 
above, let Q\ denote the cr-field generated by the random set J = {wi, . . . ,wm} and the 
values of Fx{Qw] Fj) for w ^ J . Then 



Var 



w&/Vx{Ti) 



> E Var 



^FA(g^;Fi) + ^FA(Q^;Fi 



E Var 



F\{Qw \ Tj 



using the fact that the sum over w ^ J is ^A-measurable. But by condition (ii) in 
Definition I4.lt the F\{Qw] Fj) for w G J are conditionally independent (under Qx)) so we 
obtain 



E Var 



E5^Var[FA(g^;F,)|6;,] >r7E[M], 



by condition (i) in Definition 14. 1[ Then by (14. 6p . the proof is complete. □ 



Under certain extra conditions, it is the case that 

A-iVar[r,] al (4.7) 

for some erf G [0, oo); see [4] and [11]. Often af is given explicitly as an integral; however, 
it is often non-trivial to compute or to verify that it is strictly positive. 

Under additional conditions (somewhat resembling (i) in Definition 14.11 above) it can 
be shown that af > 0. When (14.70 holds with af > for all i, we obviously have (Al). 
Conditions of this type were given in [4, 13], where a form of external stabilization is 
used (which roughly speaking says that not only do Poisson points beyond the radius 
of stabilization for x not influence x, but also x does not influence these points). The 
results of [4, 13] imply that in many cases of interest (14. 7p holds with af > (given extra 
conditions on /j and k). Functionals ^ for which this holds include those associated with 
the total edge length of the /c-nearest neighbour graph, and the total number of edges in 
the sphere of influence graph, plus others (see [4,13]). Then combining (14. 7p with external 
stabilization and the existence of moments (see Section 3 of [16] for some examples) one 
can obtain (12. 3p . 
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4.2 Example: the /c-nearest neighbour graph 

The arguments indicated above are spelled out for the particular case of the /c-nearest 
neighbour graph in Section 3.1 of [16]. Recall that for /c G M and a locally finite point 
set X C M'^, the fc-nearest neighbour (undirected) graph on X (denoted kNG(A:')) is the 
graph with vertex set A" obtained by including {x, y} as an edge whenever y & X is one 
of the k nearest neighbours of x G A", or vice versa (or both). Let ^(x; X) be one half the 
sum of the lengths in kNG(A') incident to x. Thus (for example) we have that the total 
length of kNG{X) is given by 

5^e(x;A). 

Suppose Fi, . . . , r„ are disjoint convex or polyhedral regions. We give two examples of 
conditions on {/j} and k, which, by known results together with Theorem 12.11 yield (12.31) 
for this case. 

First, suppose that k is bounded away from on UjFj. Then ^ is exponentially 
stabilizing and has moments of all orders. If fi is continuous on Fj, then (14. 7p holds with 
af > (see [16], Section 3.1). Hence Theorem 12.11 applies in this case. The conditions 
on fi and k may be relaxed (see [11]), but then extra work (such as making use of the 
nondegeneracy argument in the present paper) is needed to show that af > 0. 

Alternatively, suppose that k is equal to a positive constant on each Fj, so that 
Vx is a homogeneous Poisson point process with intensity Akj > on Fj. Suppose that 
fi = Ir^, the indicator of Fj, for each i. Then by the results of [13], we again have that 
(14. 7p holds with af > 0, and so Theorem 12. II holds. In this case, Tj is the total length of 
kNG(PAnF,). 

We conclude this section by presenting an explicit multivariate CLT of this type, 
derived from Theorem 12.11 for the case of the nearest-neighbour (directed) graph in one 
dimension. The nearest-neighbour (directed) graph on locally finite point set X is the 
graph with vertex set X obtained by including (x, y) as a (directed) edge from x E X 
to y E X when y is the nearest neighbour of x (arbitrarily breaking any ties). The 
required moments, regularity and stabilization conditions all follow from previous work 
(particularly [11,13]), and the fact that the limiting variance is non-zero follows from an 
explicit calculation (which we give below) based on the general results of [11]. 

For a finite set X C (0, 1) and a Borel set F C (0, 1), let £"(A'; F) denote the total 
weight of the nearest-neighbour (directed) graph on X, with a-power weighted edges, 
counting only edges originating from points of A" fl F. That is, if d{x; X) := d2{x] X \ {x}) 
denotes the (Euclidean) distance from x to its nearest neighbour in X , take 

i{x-X) = {d{x-X)Y, (4.8) 

for some fixed parameter a G (0, oo). Then 

For m G N, let Fi, . . . , F^ be disjoint, finite, non-null interval subsets of M. In par- 
ticular, let TTj = |Fj| G (0, oo) be the length of the interval Fj. Take fi = Ir,. Let the 
underlying density k be piecewise Borel- measurable, bounded away from and from oo, 
on each interval Fj; in particular, for each i set ) for X G Fj, where Kj G i3(Fj) 
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and Ki{x) > for all x G Fj. Consider the unmarked case (so M. = {I})- Then for A > 0, 
Vx is a Poisson point process with intensity K,i{x)X on each Fj. Using the notation of 
Theorem I2.H in this set-up we have that 

the final equality by translation-invariance. By the scaling properties ('homogeneity') of 
C, as given by (14. Sp . we have 

Ti= Y e(Aa;;APA) = A'^ J] ^(a:; Pa) = A"/:"(Pa; F,). 

All relevant stabilization, regularity and moments conditions are satisfied. Let Tii denote 
a homogeneous Poisson point process of unit intensity on (0,1), and let W„ denote a 
binomial point process consisting of n independent uniform random points on (0, 1). Then 
by Theorems 2.2 and 2.4 of [11], for a > 

lim A~Var[T,] = lim A2"-Var[/:°(pA; F,)] 

A— >oo A^oo 

= Va J Ki{x)dx + J Ki{x)dx^ , (4.9) 

where 

K := lim n2""Var[/:"(W„; (0, 1))], (4.10) 

n— >oo 

and 

5o, := E[d(0; HiT] + [ E[rf(0; Hi U {y})" - rf(0; niT]dy. (4.11) 

Let F(-) denote the (Euler) Gamma function, and let 2-Fi(', S S ■) denote the (Gauss) 
hypergeometric function (see e.g. [1], Chapter 15). By (I4.10p and equation (24) in [12], 
we have that for a > 

= (4"" + 2 ■ 3"^"^")F(1 + 2a) - 4-"(3 + a^)T{l + a)^ 

+8 ■ 6 " 'r(2 + 2a) ^p^^_^^ l + a-2 + a- 1/3). (4.12) 
(1 + a) 

We now compute 5^. By standard properties of the Poisson process, D := d{0] Hi) is 
distributed as an exponential random variable with parameter 2. So we have that for 
a > 



E[D'^] = / 2r" exp(-2r)dr = 2-"F(l + a) 
Jo 



(using Euler's Gamma integral; see e.g. 6.1.1 in [1]). By Fubini's theorem and (14. lip we 
have 



5„ = E 



D"-2 (D"-r)dt =E[D" + ((2/(l + a))-2)D^+°] 

"'0 J 

2-"F(l + a) - ^^2~i-°F(2 + a) = 2""F(1 + a){l - a), (4.13) 
1 + a 
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using the functional relation T{x) = x~^r{l+x) (see e.g. 6.1.15 in [1]) for the final equality. 
Of note is the fact that 5i = 0, so that in the a = 1 case the constant in the limiting 
(scaled) variance is the same in the Poisson and binomial cases. For a 7^ 1, 5^ > and 
the variance in the Poisson case is greater than that in the binomial case, as one expects 
(the Poisson process introduces additional randomness). 

Also by Theorem 2.1 of [11] and equation (22) in [12], we have that for a > 

X-'^E[C''{Vx;Ti)]^2-^T{l + a) [ Ki{x)dx, 

as A — >^ 00. Thus we have the following application of Theorem 12.11 

Theorem 4.1 For m G N, let Ti, . . . be disjoint intervals in M with \Ti\ = tTj G 
(0,00). Let k{x) = '^«(^)lri(a;) where, for each i, Ki G B(Ti) and Ki{x) > for all 

X G Fj. Suppose a G (0, 00). 

(i) For 1 < i < m, 



lim A"-^E[/:"(Pa; Fi)] = 2-°F(l + a) [ Ki{x)dx. 



(a) For 1 < i < m, 

lim A2°-Var[£°(PA;Fi)] = [ Ki{x)dx + ( Sa I Ki{x)dx 

where Va and Sa are given by (^T7^ and (JjJ^ respectively. 
(Hi) Given e > 0, there exists C G (0, 00) such that for all A > 1, 

..1 



sup 



P 



n 



(Var[/:°(PA;F,)])i/2 



< ti 



1=1 



<CA^-(i/2). 



Part (iii) of Theorem 14.11 is our multivariate CLT. In the particular case of piecewise 
constant k, that is Kj(x) = G (0, 00) for all x G Fj, we have that 

Ki(^X^dx = Kj|Fj| = KiTCi, 



and so, for example, af = VaUiiTi + d'^nfirf. Table [T] gives some values of the constants 
Va, given by (KT^ . and 6l, given by (KT^ . 



a 


1/2 


1 


2 


3 


4 




i + V2arcsin(l/V3) - ^ ~ 0.094148 


1 
fi 


85 


149 


135793 




7r 
.S2 





|ns 
4 


4 


4 



Table 1: Some values of V^ and 5^,. 
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